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Abstract 

We present the results of a study of the three-dimensional XY-model on 
a simple cubic lattice using the single cluster updating algorithm com- 
bined with improved estimators. We have measured the susceptibility 
and the correlation length for various couplings in the high temperature 
phase on lattices of size up to L = 112. At the transition temperature we 
studied the fourth-order cumulant and other cumulant-like quantities on 
lattices of size up to L = 64. From our numerical data we obtain for the 
critical coupling K^. = 0.45420(2), and for the static critical exponents 
7/z/ = 1.976(6) and z/ = 0.662(7). 
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I. INTRODUCTION 



Three-dimensional classical 0{N) vector models are of great interest, both as simplest sta- 
tistical models with a continuous symmetry and as a lattice version of the scalar quantum 
field theory with (0^) ^-interaction. In particular, the 3D 0(2) model, also called the XY- 
model, is relevant to the critical behaviour of a number of physical systems, such as the 
phase transition of superfluid ^He and magnetic systems with planar spin Hamiltonians. 
Quantitative knowledge of the critical behaviour of the 0{N) vector models is mostly based 
on the field theoretic renormalization group techniques at dimension D = 3 |l|] and the 
e-expansion 0,^. Very accurate values of the critical exponents are among the most suc- 
cessful predictions of these approaches. In addition, the analysis of high temperature series 
expansions |^ provides estimates for the critical temperature of particular lattice models. 

Monte Carlo simulations have also succeeded in providing detailed information about 
the critical behaviour of the 3D 0{N) vector models, but only in the case of the three 
dimensional Ising model (A^ = 1) is an accuracy close to that of analytic calculations 
reached. 

The difficulty encountered in Monte Carlo simulations with local updates is the critical 
slowing down near a phase transition. Considerable progress has been achieved during 
the last 6 years with the development of efficient non-local Monte Carlo algorithms which 
overcome critical slowing down to a large extent 0. 

In the present paper we extend previous Monte Carlo studies of the 3D XY model [^,|^ 
where cluster algorithms [llO| , pl^ were first applied to simulate this model. These studies 
were performed on vector computers, using a moderate amount of CPU-time. Since optimal 
vectorization cannot be reached for cluster algorithms, we used modern RISC stations for the 
present study. Using about two months of CPU-time we were able to simulate larger lattices 
and reached a considerably better statistical accuracy than in the previous studies. This 
accuracy allowed us to control systematic errors in the estimates for the critical coupling 
and the critical exponents. 

This paper is organized as follows. In section 2 we give the definition of the model 
and describe the cluster updating algorithm, in section 3 we give our results for the high 
temperature phase, while section 4 contains our data obtained in the critical region. In 
section 5 we compare our results with those of previous studies. 



II. CLUSTER UPDATE MONTE CARLO OF THE SD XF-MODEL 

We study the XF-model in three dimensions defined by the partition function 

Z = Y[ dsi exp(ir ^ Si ■ Sj) , (1) 

where Sj is a two dimensional unit vector, the summation is taken over all nearest neighbour 
pairs of sites i and j on a simple cubic lattice A and K = is the coupling, or more 
precisely, the reduced inverse temperature. 
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For ferromagnetic interactions J > 0, the Xy-model has a second-order phase transition 
separating a low temperature phase with non-zero magnetization from a massive disordered 
phase at high temperatures. This phase transition can be viewed alternatively as due to 
Bose condensation of spin waves [|12| or the unbinding of vortex strings [0,0. 



A major difficulty encountered in Monte Carlo simulation at second-order phase transi- 
tions is critical slowing down. The autocorrelation time r, which is roughly the time needed 
to generate statistical independent configurations, grows as r oc at criticality, where L 
is the linear size of the system and z is the dynamical critical expontent. Random walk 
arguments indicate that local updates like the Metropolis algorithm result in z = 2, which 
is consistent with the numerical finding for the 3D XF-model 

In the case of 0{N) vector models, critical slowing down can be drastically reduced. 



using cluster algorithms [H,ll 



In the present work we employ the single cluster algorithm which was introduced by 



Wolff |Tl|. Let us shortly recall the steps of the update. First choose randomly a reflection 
axis in the plane. Denote the component of the spin Si which is parallel to this reflection 
axis by sf and that which is orthogonal by s^. Then choose randomly a site i of the lattice 
as a starting point for the cluster C. Visit all neighbour sites j of i. These sites join the 
cluster with the probability 

j) = 1 - exp(-i^(.^^^ + \sisf\)). (2) 

After this is done, visit the neighbours of the new sites in the cluster and add them to the 
cluster with probability p{i,j) which is given above. Iterate this step until no new sites enter 
the cluster. Now flip the sign of all S"*" contained in the cluster. 



III. NUMERICAL RESULTS IN THE HIGH TEMPERATURE PHASE. 

A. Observables to be measured 

Let us flrst summarize the definitions of the observables that we studied. The energy 
density is given by the two-point correlation function G{xi,Xj) = (siSj) at distance one 

E = ^,T.{^^s,). (3) 

The specific heat of the system at constant external field is defined by the derivative of 
the energy density with respect to the inverse temperature. It can be obtained from the 
fiuctuations of the energy H = — SiSj 

Ch = ^ {{H') - (Hf) . (4) 

The magnetic susceptibility x gives the reaction of the magnetization m = J2ieA to 
external field. In the high temperature phase one gets 
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(5) 



since (m) = 0. 

Cluster algorithms enable one to reduce the variance of the expectation values in the 



high temperature phase by using improved estimators [^,0. The improved estimator of 
the magnetic susceptibility is given by 



X.VI 



mp 



i&C 



(6) 



where \C\ denotes the number of spins in the cluster C. 

There are two common definitions of a correlation length ^. The exponential correlation 
length ^exp is defined via the decay of the two-point correlation function at large distances 



exp 



lim 



log G{xi,Xj 



(7) 



which is equal to the inverse mass gap. For the measurement of the exponential correlation 
length we consider the correlation function 



G(t) = {OoOt) oc fexp(^)+exp( 



exp 



(8) 



exp 



of the translational invariant time slice magnetization Ot = J2i s{xi,t). 
The second-moment correlation length is defined by 



^2 



(Mil 



V4sin2(7r/L) 



(9) 



with F = G{k)\\k]=2TT/L, where G{k) = I]jgA(exp(iA;Xj)soSj) is the Fourier transform of the 
two-point correlation function and x the magnetic susceptibility. For more details see for 
example ref. |l7|. The two definitions of the correlation length do not coincide, since in ^^xp 
only the first excited state enters, while in the case of .^2"'* a mixture of the full spectrum is 
taken into account. However, near the critical point the two quantities should scale in the 
same way. As for the magnetic susceptibility there exist improved estimators for the two 
definitions of the correlation length. The improved estimator of the two-point correlation 
function is given by 



imp 



(10) 



where Sij{C) = 1 if z and j belong to the same cluster C, otherwise Sij{C) = |[15|J16[| . For ^2"^ 
one has to provide a Fimp- This is given by the Fourier transform of the improved two-point 



correlation function 



imp 



G{k) 



imp 



2n ' 
" L 



(11) 
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with 

The hehcity modulus describes the reaction of the system to a suitable phase twisting field 




Y: si cos{kx,)] + si fAnikxA I \ . (12) 

i&C I \i€C /I 



T9| . The lattice definition of the helicity modulus is given by 



where /i is a unit vector in x,y or z direction and e^jj) the unit vector connecting the sites 



i and j 



B. Monte Carlo Simulations 

In order to obtain an estimate of the critical coupling and determine static critical 
exponents, we have done 15 simulations at couplings from K = 0.4 up to = 0.452 on 
lattices of linear size L = 24 up to L = 112. The simulation parameter and the results of the 
runs are given in Tables | and ||. The statistics is given in terms of measurements taken 
every Nq update steps. A^o is chosen such that approximately Nq x (C) = and hence the 
whole lattice is updated once for a measurement. We estimated the statistical errors aA of 
expectation values {A) from 

= N/{2r) 

and from a binning analysis. These error estimates were consistent throughout. The statisti- 
cal error of quantities which contain several expectation values we calculated from Jackknife- 



blo eking [21 . 



1. Finite-size effects 

We tried to avoid a finite-size scaling analysis. Hence we had to choose our lattices large 
enough to ensure that deviations of the values of the observables from the thermodynamic 
limit values are negligible. 

We therefore have measured the energy density the specific heat C/i, the helicity 
modulus T , the exponential correlation length ^exp and the second-moment correlation 
length ^2"d for fixed coupling K = 0.435 and increasing system size L = 4 up to L = 32. 



The results are summarized in Table |n|. The values of the observables obtained for L = 24 
and L = 32 are consistent within error bars. Furthermore, the values of the helicity modulus 
T for L = 24 and L = 32 are consistent with 0, which is the thermodynamic limit value of 
the helicity modulus in the high temperature phase. The correlation length at = 0.435 
is approximately 4. Hence we conclude, assuming scaling, that the systematical deviations 
from the thermodynamic limit are smaller than our statistical errors for L/^ > 6. This 
condition is fulfilled by all the simulation parameters of our runs given in Table |. 
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2. Energy density and specific heat 

The energy density E shows, as expected, no singular behaviour close to the critical 
temperature. In the scaling region the specific heat Ch should follow 

Ch = Creg + Co j ) (15) 

where Creg denotes the regular part of the specific heat and a is the critical exponent of the 
specific heat. In order to estimate a we did a four-parameter least-square fit. However, it was 
not possible to extract meaningful estimates. The best fit to the data leads to Kc = 0.456(2), 
a = 0.23(13) with x^/d.o.f. ~ 0.93 and the relative errors of the constants are about 100%. 
If we fix the critical coupling to Kc = 0.45420 (this is our estimate obtained at criticality) 
the quahty of the fit gets worse. Therefore we assumed a = and fitted the data following 

Ch = Creg + Cq log (^^^^) • (16) 

The best three-parameter fit to our data leads to Kc = 0.4543(4), Creg = —0.49(20) and 
Cq = —1.61(7) with x^/d.o.f. ~ 0.61, where data with ^ > 2.5 are taken into account. 
This result shows that our data for the specific heat, combined with extended ansatze, are 
compatible with an a = —0.007(6) obtained from the hyperscaling relation a = 2 — Du and 
the estimate v = 0.669(2) from resummed perturbation series |T[], but have no predictive 
power for the exponent a. 



3. Magnetic susceptibility 

For comparison we give in Table | the results for the standard and the improved sus- 
ceptibilities. The statistical error of Ximp is about 3.5 to 8 times smaller than the error of 
the standard susceptibility. But one should remark that the statistical error of the standard 
estimator depends very much on how often one measures. In the following we only discuss 
the results obtained with the improved estimator. In order to estimate the critical coupling 
and the susceptibility exponent 7 we performed a three-parameter least-square fit following 
the scaling law 

f Kc-K y^ 

X = Xo — • (17) 



K 



c 



We obtained 7 = 1.324(1), Kc = 0.454170(7) and Xo = 1-009(2) with x^/d.o.f. = 0.65, when 
all data are taken into account. In order to test the stability of the results we successively 
discarded data points with small K. The results of these fits are summarized in Table 
X^/d.o.f. remains small and the results for 7, Kc and Xo are consistent within the error 
bars for all data-sets that we used. But the small x^ /d.o.f. of the fits discussed above, is of 
course, no proof for the absence of corrections to the scaling. From renormalization group 
considerations |2^ one expects confiuent and analytical corrections of the type 
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XiK) = XO ^ j + Xconf. ^ j + Xanal. ^ j , (18) 

with Ai = uju, where u is the critical exponent of the correlation length and uj denotes 
the correction-to-scaling exponent. We fitted our data according to the scaling law with 
corrections. Since a fit with 6 free parameters is hard to stabilize, we fixed the critical 
exponents to the values 7 = 1.3160(25), u = 0.780(25) and z/ = 0.669(2) which are obtained 
from resummed perturbation series [|I|. Including all the data points in the fit we get Kc = 
0.454162(9), xo = 1-058(7), Xconf. = -0.16 = (6) and Xanai. = 0.18(10) with xVd.o.f. ^ 
0.73. The Xo which is obtained from the simple scaling fit (|T7|), and that obtained from the 
fit allowing corrections to the scaling, differ by a larger amount than their statistical errors. 
This shows that one cannot interpret a small x^/d-o.f. as the absence of systematic errors 
due to an incomplete fit ansatz. 

One can also write the scaling relations in terms of the temperature T = ^. This leads 

to 

m = xo[^^) (19) 



X{T) = Xo + Xconf. + Xanal. (20) 



and 



rj-i J • /\,t^uiLj . r-j-\ J ' /\.UiiLu,L. rj-i 

with corrections. We repeated the analysis as done above. Taking all 15 data points into 
account we get for the simple scaling fit a x^ jd.o.f ~ 61.2. We again subsequently discarded 
data points with small K. A summary of the results is given in Table 0. Starting from 5 
discarded data points the x^jd.o.f. becomes approximately 1. But the results obtained for 
Kf-i 7 and Xo ^^rc not consistent with those obtained from the fit according to eq. (p!^). 

Finally, we performed a four-parameter fit to the scaling relation with corrections and 
fixed values for the exponents. Taking all data points into account we obtain Kc = 
0.454163(9), Xo = 1-059(7), Xconf. = -0.17(6) and Xanai. = 1-59(10) with xVd.o.f. ^ 0.71. 
The results for Xo and Xconf of the fits according to the ansatze (p!8| ) and (^) are consis- 
tent within the error bars. The ambiguity between the ansatz with K as variable and that 
with T as variable is covered by the analytic corrections. 

We conclude that the scaling ansatz (|1^ fits well if one chooses the coupling K as the 
variable. But we also learned that a small x^ /d.o.f. does not exclude systematic errors, due 
to corrections to the scaling, that are larger than the statistical ones. Hence it is hard to give 
final estimates obtained from the simple scaling ansatz that also take systematic errors into 
account. From the ansatz with corrections to the scaling we obtain, assuming 7 = 1.3160(25) 
and Ai = 0.52182, the results = 0.454162(13), Xo = 1-058(22) and Xconf. = -0.16(11), 
where the uncertainty of 7 is taken into account. 

We also like to emphasize that the Wegner amplitude Xconf. is negative for the fits that 
take corrections into account. This is in agreement with a field-theoretical renormalization 
group calculation of Esser and Dohm, which predicts the confluent correction-to-scaling 
amplitude to be negative for a finite cut-off ||23l . 
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^. Correlation length 



We extracted ^exp from the large distance behaviour of the improved time shce correlation 
function eq.(^. We therefore considered the effective correlation length, defined by 

where for brevity we have suppressed the contribution due to periodic boundary conditions. 
As an example, we show in Fig. |l| the results for ^effif) obtained on a 1123 lattice at 
K = 0.452. A single state dominates the correlation function and a plateau sets in around 
t = and extends to t = S^exp, with no visible degradation due to increasing statistical 

errors at large t. As our final estimate for ^^xp we took self-consistently ^eff(t) at the distance 

In order to calculate ^2"-'' we used the improved version of eq.@. The advantage of this 
definition is that no fit is needed to obtain the correlation length. The data of C,exp and ^2"- 



are given in Table ||. The deviation of ^2"'* from C,exp is about 1% for K = 0.40 and becomes 
smaller than 0.1% for > 0.448. 

Since the difference of ^exp and ^2"'* is so small, we will discuss only the results of ^exp in 
the following. First we did a three-parameter fit for C,exp following the simple scaling ansatz 

/K - K\ 



The results are given in Table Vl. Taking all data into account, the fit has a large y^/d.o.f 
of about 9. Starting with three data points with small K being discarded, the x^ld-o.f. is 
close to 1. But still the critical exponent v and the critical coupling systematically tend to 
smaller values. 

If we fit the data to the simple scaling ansatz (^2]), where the coupling is replaced as 
variable by the temperature, a similar behaviour is observable. The results are shown in 
Table |VII| . Here one also has to discard three data points to obtain a /d-o.f. close to 1. 
But now the estimates of the critical exponent and the critical coupling start at lower values 
and tend to larger ones. 

This indicates that corrections to the simple scaling ansatz have to be taken into account. 
Therefore we have fitted all data to the scaling relation with corrections given by 

{'^r\ (.3) 

whereas, in the case of the magnetic susceptibility, Ai = ujv. The four-parameter fit to all 
data points with the critical exponents fixed to the resummed perturbation series estimates 
leads to = 0.454167(10), = 0.498(2), ^^onf. = -0.10(2) and ^anai. = -0.07(4) with 
X^/d.o.f. ^ 0.63. 

We also made a four-parameter fit to the scaling relation with corrections where the 
coupling is replaced by the temperature. This leads to Kc = 0.454165(10), = 0.498(2), 
lonf. = -0.09(2), ianai. = 0.24(3) with xVd.o.f. ^ 0.63. 
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Taking the uncertainty of v into account, we leave at Kc = 0.454166(15), = 0.498(8) 
and ^conf. = -0.10(6). 

In summary, we conclude that systematic deviations from the simple scaling ansatz ( p2D 
due to corrections to scaling are important for the analysis of the correlation length data 
in the coupling range that is accessible to Monte Carlo simulations. Thus it is hard to 
obtain accurate estimates for the critical exponents and the critical coupling from such an 
approach. 



IV. NUMERICAL RESULTS AT CRITICALITY 



On lattices of the size L = 4, 8, 16, 32 and 64 we performed simulations at Kq = 0.45417 
which is the estimate for the critical coupling obtained in the previous section. As in the 
high temperature simulation the statistics are given in terms of N measurements taken 
every A^o update steps. We have chosen Nq such that on the average the lattice is updated 
approximately twice for one measurement. The results of the runs are summarized in Table 



A. Phenomenological Renormalization Group 



First we determined the critical coupling and the critical exponent ly employing 
Binder's phenomenological renormalization group method p4[. In addition to the fourth- 
order cumulant defined on the whole lattice we also studied cumulants defined on subblocks 
of the lattice. Therefore let us first introduce blockspins 

5B = Lr"-^^ AE^1> (24) 

where is the size of the block and 1/2{D — 2) is the canonical dimension of the field. In 
particular we studied the fourth-order cumulant 



< (Sir > 
'I 

for Lb = L,L/2 and a nearest neighbour interaction on subblocks 

< SbiSb2 > 



NN = ^-y^ (26) 



for Lb = L/2. 

For the extrapolation to couplings K other than the simulation coupling Kq, we used 
the reweighting formula 

E.^^exp((-ir + Ko)g,) 
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where i labels the configurations generated according to the Boltzmann- weight at Kq. We 
computed the statistical errors from Jackknife binning on the final result of the extrapolated 
cumulants. The extrapolation only gives good results within a small neighbourhood of the 
simulation coupling Kq. This range shrinks with increasing volume of the lattice. 

For sufficiently large Lb the cumulants have a non-trival fixed point at the critical cou- 
pling. When one considers the cumulants as a function of the coupling, the crossings of the 
curves for different L provide an estimate for the critical coupling K^. As an example we 
show in Fig. ^ the fourth-order cumulant in a neighbourhood of Kc- The figure shows that 
the crossings of the cumulant are well covered by the extrapolation (|27|) . The error bars of 
Ul with L = 64 blow up for \K - Ko\ > 0.001 while \K^ross - ^o| = 0.00003 for L = 32 
and L = 64. The results for the crossings are summarized in Table The given errors 
are taken from the size of the crossings of the error bars. The convergence of the crossing 
coupling Kcross towards Kc should follow 

KcrossiL) = Kc{l + const.L-^^+^^), (28) 



where u is the correction to scaling exponent |2^. Our data for the crossings of the cumulants 
did not allow us to perform a two-parameter fit, keeping the exponents fixed, following the 
above formula. Within the statistical errors the results of the crossings of the fourth-order 
cumulants on L = 8 and L = 16 up to L = 32 and L = 64 are consistent. The convergence of 
the crossings of NN towards Kc seems to be slower than that of the fourth-order cumulants, 
but it is interesting to note that the Kcross for the fourth-order cumulant and NN come 
from different sides with increasing L. This is shown in Fig. ^ where the estimates of 
Kcross versus the lattice size L are plotted. Our final estimate for the critical coupling is 
Kc = 0.45420(2) obtained from the L = 32 and L = 64 crossing of the fourth-order cumulant 
on the full lattice. Taking into account the fast convergence of the crossings towards Kc, 
that is predicted by (|28D , we conclude that the systematic error of our estimate for Kc is 
smaller than the given statistical error. 

At the critical coupling Kc the cumulants converge with increasing lattice size L to a 
universal fixed point. The convergence rate is given by 

ULiKc) = U^il + const.L'^) . (29) 

The results for the cumulants aX K = 0.45420 , which is our estimate of critical coupling, 
are given in Table The data did not allow us to perform a two parameter fit with uj being 
fixed. Hence we take the value Ul{Kc) = 0.589(2) form L = 64 as our final estimate for the 
fixed point of the fourth-order cumulant on the full lattice, where we now have taken into 
account the uncertainty of the estimated critical coupling. 

We extracted the critical exponent u of the correlation length from the L dependence of 
the slope of the fourth-order cumulant at criticality According to Binder, the scaling 



relation for the slope of the fourth-order cumulant is given by 

dU{L,K) 



dK 



oc L'/". (30) 
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We evaluated the slopes of the observables A entering the cumulant U according to 

d{A) 



dK 



{AH)-{A){H), (31) 



where A is an observable and H is the energy. The statistical errors are calculated from 
a Jackknife analysis for the value of the slope. First we estimated the exponent u from 
different lattices via 

. . , .■"<^;)-'°<^;'.., (32) 



In 



fdA{L2,K) 



fdAiL,,K) 



The results are given in Table The estimates for v stemming from Ul and f/L/2 are stable 
with increasing L and consistent with each other for L > 16. Therefore we performed a fit 
according to eg. (pO|) with Ul from lattices of the size L = 16 up to L = 64. We consider the 



result u = 0.662(7) as our final estimate for the critical exponent of the correlation length. 

B. Magnetic Susceptibility 

In order to estimate the ratio 7/z/ of the critical exponents we studied the scaling be- 
haviour of the magnetic susceptibility defined on the full lattice and on subblocks. The 
dependence of the susceptibility on the lattice size is given by 

X « L-^^" (33) 

at the critical coupling. We have estimated 7/z/ from pairs of lattices with size Li, L2. The 
ratio is then given by 

V In(Li) - ln(L2) ' ^ 



The second column of Table PCII| shows the estimates of the ratio obtained from the sus- 
ceptibility defined on the full lattice, while the third column shows the estimates obtained 
from the blockspin-susceptibility with subblocks of the size L/2. The estimates for 7/z/ ob- 
tained from the subblocks monotonically increase with increasing lattice size L, while those 
obtained from the full lattice decrease. The results obtained from the full lattice for L > 16 
and the result from the subblocks of the largest lattices are consistent within error bars. 
Hence we take 7/z/ = 1.976(6) as our final result, where statistical as well as systematic 
errors should be covered. Using the scaling relation 17 = 2 — ^ we obtain for the anomalous 
dimension r] = 0.024(6). 

C. Hyper scaling and specific heat 



In ref. p4[ a dangerous irrelevant scaling field u is proposed as explanation for a possible 



violation of hyperscaling. Dangerous means that the scaling function of the correlation 
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length vanishes with some power q of the vanishing irrelevant scaling field. Hence the 
correlation length should scale as 



i oc L^+'^y- (35) 

at the critical point. Remember that Uu is negative for an irrelevant scaling field. At 
Kc = 0.45420, which we obtained from the analysis of the fourth-order cumulant, we have 
fitted ^2"rf to this relation. The reweighted estimates of ^2"^ are shown in Table PCII1| . Taking 



lattices of size L = 16 up to L = 64 into account we estimate qi/u = 0.007(2) with x^/d.o.f. ~ 
0.3 and only statistical errors considered. This indicates that there is no or only very small 
hyperscaling violation due to a dangerous irrelevant field. 
At criticality the specific heat should scale as 

Ch{L) = Creg + COUSt. L~ , (36) 

where Creg denotes the regular part of the specific heat and a is the critical exponent of 
the specific heat. Using the critical exponent u = 0.662(7) obtained from the analysis above 
the hyperscaling relation a = 2 — Dv gives a = 0.014(21). We also tried to estimate a via 
a three-parameter fit, following the finite-size scaling relation. However, we are not able to 
give a stable estimate for a. 

D. Helicity modulus 

The 3-D XY model is assumed to share the same universality class as an interacting Bose 
fluid, and the helicity modulus T should be proportional to the superfluid density Qs of the 
Bose fluid [0. Near the critical coupling the superfluid density, resp. the helicity modulus 
should scale as 

Qs (X T (X \K — Kc\" , (37) 

with V the critical exponent of the superfluid density. Assuming hyperscaling the Josephson 
relation reads v = {D — 2)v [|l^. Hence the product 



T ■ L = const (38) 

should stay constant at the critical point in 3D. To check this prediction we have measured 
the helicity modulus T on lattices of size L = 4 to L = 32. The estimator of TL becomes 
noisy with increasing lattice size. We tried to overcome this problem by measuring more 
often, which did not remove the problem completely. Hence we skipped the measurement of 
T for L = 64. The results, shown in Table pC114 indicate that the above relation holds. 

E. Performance of the Algorithm 

The efficiency of a stochastic algorithm is characterized by the autocorrelation time 
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-| oo 

r = ^ E Pit) ^ (39) 
where the normahzed autocorrelation function p{t) of an observable A is given by 

n(t) = (-4^ • Aj+t) - {Af 

^ ' (v42) - {Af ■ ^ ' 

We calculated the integrated autocorrelation times r with a self-consistent truncation win- 
dow of width 6r for the energy density E and the magnetic susceptibility x fo^' lattices with 
L = 4 up to L = 64 at the coupling K = 0.45417. In Fig. ^ we show a log-log plot of the 
integrated autocorrelation times r of the energy density E and the magnetic susceptibility 
X versus the lattice size L given in units of the average number of clusters that is needed 
to cover the volume of the lattice. Our estimates for the critical dynamical exponents are 
ze = 0.21(1) and = 0.07(1) taking only statistical errors into account. These results are 
consistent with those of Janke 0. 

Finally let us briefly comment on the CPU time: 160 single cluster updates of the 64^ 
lattice at the coupling K = 0.45417 plus one measurement of the observables took on average 
26 sec CPU time on a IBM RISC 6000-550 workstation. All our MC simulations of the 3D 
XY model together took about two months of CPU-time. 



V. COMPARISON OF OUR RESULTS WITH PREVIOUS STUDIES 



In Table pClV| we display estimates of critical properties of the 3D XF-model obtained by 
various methods. Our estimates of from the scaling fit to the high temperature data and 
from the penomenological RG approach are consistent within 2 standard deviations. But 
only for the result from the phenomenological RG approach are the systematical errors fully 
under control. Our error of Kc is about 4 times smaller than that of previous MC studies 
PJ^, and also about 4 times smaller than that obtained recently p6[ from the analysis of 
a 14:th order high temperature series expansion [28|. Recently Butera et al. extended 



the high temperature series expansion for the sc lattice to the order 17. Their value for the 
critical coupling is by three times their error estimate smaller than our value. 

The error of 7 obtained from a fit to the simple scaling ansatz is about 5 times smaller 
than those of previous MC studies Pfl; however, the systematical errors are not under 
control. The value of 7 is, within two standard deviations, consistent with the estimate of 
Ref . . Our estimate of 7 is consistent with the value obtained from the high temperature 
EB|,|2B|,B7| and, within two standard deviations, consistent with the value of 



series expansion 

the e-expansion while the very accurate estimate from the resummed perturbation series 
|jl|] is smaller than our estimate by three times the quoted error. 

Our estimate for u is consistent within error bars with all other estimates we quote in 
Table pCIV| . Our quoted error bars are 3.5 times larger than that of ref. @. Janke used finite 



differences to determine the slope of the cumulant [29|, while we used fluctuations at a single 
temperature (|3lD . Furthermore the smallest lattice size L = 16 included in our fit is chosen 
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to be rather conservative. The most accurate number for v stems from the measurement of 
the superfluid fraction of "^He |^ . 



In this work we give for the first time an accurate direct MC estimate for the exponent 
7], the anomalous dimension of the field. The uncertainty of the estimate is comparable with 
those obtained with field theoretical methods. Our value of the exponent r] is consistent 
with the estimates from the high temperature series expansion p5| , |26| and with that of the 
resummed perturbation series []I| , but is smaller than the e-expansion result by more than 
twice our error estimate. 

Our result for the critical fourth-order cumulant, is consistent with previous MC results 
PJ^. But the value obtained from e-expansion |^ is off by about 20 times our error estimate 
that also takes into account systematic errors. Furthermore we provide estimates for the 
critical fourth order cumulant on subblocks and a nearest neighbour blockspin product AW. 
These numbers might be useful in testing other models sharing the XY universality class. 



VI. CONCLUSIONS 



The application of the single cluster algorithm which is almost free of critical slowing 
down for the 3D XY model, and the extensive use of modern RISC workstations allowed 
us to increase the statistics as well as the studied lattices sizes considerably compared with 
previous MC simulations In the high temperature phase of the model we measured 

correlation length up to 17.58 with an accuracy of about 0.1%. But the analysis of our data 
for the correlation length and the magnetic susceptibility showed that it is hard to control 
systematic errors due to confluent and analytic corrections. It seems to be much easier 
to fight the systematic errors in the phenomenological RG approach. Analytic corrections 
are absent at the critical point and corrections to the scaling are less harmful, since the 
relevant length scale at criticality is the lattice size, which can be chosen much larger than 
the correlation length in the thermodynamic limit of the high temperature phase. From the 
crossings of the fourth-order cumulant we obtain = 0.45420(2), which reduces the error 
by a factor of about 4 compared with previous MC studies PJ^. Further improvements of the 
accuracy of the estimates of the critical coupling and critical exponents seem to be reachable 
by just increasing the statistics, while keeping the present lattice sizes. The accurate values 
obtained for critical cumulants could be very useful for testing whether other models share 
the XY universality class. Here of course a proper block-spin definition is essential. 
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FIGURES 



FIG. 1. The effective correlation length Ce//(i) as a function of separation for couphng 
K = 0.452 on a lattice of size L = 112. 

FIG. 2. Reweighting plot of the Binder cumulant Ul oi the full lattice from the simulation at 
K = 0.45417. The dashed lines give the statistical errors obtained by a binning analysis. 

FIG. 3. Plot of the convergence of the critical coupling obtained by the cumulant crossing 
method. Because of the small statistical errors one is able to see systematic convergence of the 
critical coupling. 

FIG. 4. Intergated autocorrelation times r of the energy density E and the magnetic suscep- 
tibility X versus the lattice size L. The dynamical critical exponent is given by the slopes of the 
fits. 
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TABLES 



TABLE I. Results of the energy density E, the specific heat at constant external field C^, the 

impovcd susceptibility Ximp and the standard susceptibility x obtained from the simulations in the 
high temperature phase. The parameters of the runs are given in terms of the simulation coupling 
A, the linear size of the system L, and the statistics, with N the number of measurements taken 
every Nq update steps. 



K 


L 


N 


Ao 


E 


Ch 


Ximp 


X 


0.400 


24 


20k 


0.9k 


0.24697(6) 


3.17(3) 


16.848(14) 


16.70(12) 


0.410 


32 


10k 


2k 


0.25779(5) 


3.39(6) 


22.108(19) 


22.24(32) 


0.420 


32 


10k 


1.6k 


0.26970(5) 


3.65(6) 


30.994(29) 


31.18(32) 


0.425 


32 


10k 


1.2k 


0.27605(6) 


3.96(7) 


38.21(5) 


38.47(38) 


0.430 


32 


20k 


0.8k 


0.28286(4) 


4.22(4) 


49.12(6) 


48.87(36) 


0.435 


32 


20k 


0.6k 


0.29023(7) 


4.57(8) 


66.57(15) 


66.6(7) 


0.437 


32 


20k 


0.4k 


0.29329(5) 


4.72(6) 


77.11(16) 


76.9(5) 


0.440 


32 


20k 


0.4k 


0.29818(5) 


5.11(6) 


99.50(20) 


99.2(7) 


0.443 


48 


21k 


Ik 


0.30336(3) 


5.39(7) 


136.41(20) 


135.8(9) 


0.445 


48 


40k 


0.5k 


0.30705(3) 


5.84(6) 


177.04(30) 


176.2(9) 


0.448 


64 


25k 


Ik 


0.31303(2) 


6.43(8) 


299.14(53) 


301.0(2.0) 


0.449 


64 


25k 


Ik 


0.315214(18) 


6.61(7) 


377.75(64) 


378.9(2.5) 


0.450 


96 


12k 


3k 


0.317497(15) 


7.02(11) 


503.85(85) 


506.5(4.6) 


0.451 


96 


12k 


2k 


0.319854(13) 


7.33(10) 


722.1(1.3) 


724.6(6.5) 


0.452 


112 


12k 


2k 


0.322446(12) 


8.10(14) 


1193.0(3.0) 


1201.0(10.0) 
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TABLE II. Results for the correlation lengths ^exp and ^2"'* obtained from the simulations in 
the high temperature phase. i2nd/£,exp gives the ratio between the values of the two correlation 
lengths. 



K 




L 






^exp 




^2"'* /Cexp 


0.400 




24 




1.876(1) 


1.898(2) 




0.9884 


0.410 




32 




2.182(1) 


2.202(2) 




0.9909 


0.420 




32 




2.624(2) 


2.639(2) 




0.9943 


0.425 




32 




2.938(2) 


2.953(2) 




0.9949 


0.430 




32 




3.361(2) 


3.375(3) 




0.9959 


0.435 




32 




3.947(5) 


3.959(5) 




0.9969 


0.437 




32 




4.262(5) 


4.273(6) 




0.9974 


n /I /in 

U.44U 








4.875(5) 


4.000(^0 j 




n QQ7Q 


0.443 




48 




5.746(5) 


5.756(6) 




0.9983 


0.445 




48 




6.582(6) 


6.593(7) 




0.9983 


0.448 




64 




8.638(9) 


8.645(10) 




0.9991 


0.449 




64 




9.738(9) 


9.747(10) 




0.9991 


0.450 




96 


11.288(10) 


11.295(10) 




0.9994 


0.451 




96 


13.587(16) 


13.594(16) 




0.9995 


0.452 




112 


17.570(19) 


17.580(20) 




0.9994 


TABLE III. Results of the energy density the specific heat C/j, the helicity modulus T, the 


improved susceptibility Ximp^ the correlation lengths ^(.^p and .^2"'* a,nd the ratio L/^2"'* 


obtained 


from simulations 


at K = 0.435 with linear system size 


L. The statistics is j 


fiven in terms of N 


measurements taken every Nq update steps. 








L N 


No 


E 


Ch 


T 


Ximp ^2^<^ 


^exp 




4 25k 


0.1k 


0.3621(26) 


6.43(13) 


0.2180(44) 


15.83(21) 2.015(19) 


2.08(2) 


1.98 


8 20k 


0.2k 


0.3054(26) 


6.56(13) 


0.0608(16) 


43.00(07) 3.228(3) 


3.272(4) 


2.48 


16 20k 


0.2k 


0.29114(9) 


4.92(6) 


0.0056(17) 


64.31(16) 3.886(6) 


3.901(6) 


4.12 


20 20k 


0.3k 


0.29049(7) 


4.71(5) 


0.0031(16) 


66.21(14) 3.938(4) 


3.949(5) 


5.08 


24 20k 


0.3k 


0.29024(7) 


4.65(5) 


-0.0001(17) 


66.55(15) 3.946(6) 


3.959(6) 


6.08 


32 10k 


0.6k 


0.29013(7) 


4.52(9) 


0.0009(25) 


66.64(15) 3.949(5) 


3.960(5) 


8.10 
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TABLE IV. Estimates for the critical coupling Kc, the static critical exponent 7 and the 
amplitude xo obtained from a fit of the improved susceptibility Ximp to eg. (p^ . x^/d-o.f gives the 
quality of the fit. # denotes the number of discarded data points at small couplings. 



# 


Kc 


7 


Xo 


xVd. 





0.454170(7) 


1.3241(10) 


1.0090(24) 


0.65 


1 


0.454168(8) 


1.3238(12) 


1.0099(30) 


0.69 


2 


0.454175(9) 


1.3252(14) 


1.0057(40) 


0.49 


3 


0.454173(10) 


1.3248(18) 


1.0069(53) 


0.53 


4 


0.454170(11) 


1.3239(22) 


1.0101(66) 


0.52 


5 


0.454179(14) 


1.3264(31) 


1.0016(100) 


0.41 


6 


0.454176(15) 


1.3256(35) 


1.0043(117) 


0.44 


7 


0.454174(17) 


1.3251(43) 


1.0060(145) 


0.52 


8 


0.454174(19) 


1.3250(52) 


1.0063(180) 


0.65 


9 


0.454180(24) 


1.3272(74) 


0.9979(267) 


0.81 


10 


0.454197(42) 


1.3337(156) 


0.9736(576) 


1.11 


11 


0.454208(58) 


1.3384(240) 


0.9557(885) 


2.14 



TABLE V. Estimates for the critical coupling K^, the static critical exponent 7 and the am- 
plitude Xo obtained from a fit of the improved susceptibility to eg. (|l9|) . x^ jd.o.j gives the guality 
of the fit. denotes the number of discarded data points at small couplings. 



# 


Kc 


7 


Xo 


X^ld.o.f. 





0.453871(6) 


1.2351(8) 


1.3972(27) 


77.1 


1 


0.453932(7) 


1.2471(10) 


1.3500(34) 


40.0 


2 


0.453995(8) 


1.2604(13) 


1.2970(45) 


12.8 


3 


0.454028(9) 


1.2683(16) 


1.2648(59) 


6.98 


4 


0.454040(10) 


1.2733(20) 


1.2440(74) 


5.27 


5 


0.454091(13) 


1.2849(29) 


1.1959(111) 


1.39 


6 


0.454098(14) 


1.2871(33) 


1.1870(128) 


1.31 


7 


0.454110(16) 


1.2907(40) 


1.1720(158) 


1.08 


8 


0.454120(18) 


1.2938(49) 


1.1587(196) 


1.02 


9 


0.454135(23) 


1.2991(70) 


1.1359(289) 


0.99 


10 


0.454166(42) 


1.3112(150) 


1.0845(618) 


1.06 


11 


0.454182(59) 


1.3177(231) 


1.0570(945) 


1.98 
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TABLE VI. Estimates for the critical coupling Kc, the static critical exponent v and the 
amplitude .^o obtained from a fit of the exponential correlation length ^exp to eq.(p2|). /d.o.f 
gives the quality of the fit. 7^ denotes the number of discarded data points at small couplings. 



# 


Kc 


V 


eo 







0.454325(9) 


0.7029(7) 


0.4294(8) 


9.2 


1 


0.454301(9) 


0.7003(8) 


0.4327(9) 


4.9 


2 


0.454286(10) 


0.6985(9) 


0.4351(11) 


3.6 


3 


0.454269(11) 


0.6964(10) 


0.4381(13) 


2.0 


4 


0.454247(12) 


0.6933(13) 


0.4426(18) 


0.66 


5 


0.454243(14) 


0.6927(18) 


0.4436(25) 


0.71 


6 


0.454235(16) 


0.6914(21) 


0.4457(31) 


0.59 


7 


0.454223(18) 


0.6895(24) 


0.4487(37) 


0.31 


8 


0.454216(20) 


0.6882(30) 


0.4509(48) 


0.25 


9 


0.454210(24) 


0.6870(39) 


0.4529(64) 


0.26 


10 


0.454208(46) 


0.6866(90) 


0.4537(157) 


0.39 


11 


0.454218(66) 


0.6890(140) 


0.4493(246) 


0.74 



TABLE VII. Estimates for the critical coupling K^., the static critical exponent v and the 
amplitude ^0 obtained from a fit of the exponential correlation length ^^^.^p to eq. (|2^) , where the 
coupling is replaced by the inverse temperature. jd.o.f gives the quality of the fit. # denotes 
the number of discarded data points at small couplings. 



# 




V 


Co 







0.454079(8) 


0.6612(6) 


0.5024(8) 


8.5 


1 


0.454098(8) 


0.6632(7) 


0.4994(9) 


5.3 


2 


0.454118(9) 


0.6656(8) 


0.4957(11) 


1.7 


3 


0.454125(10) 


0.6665(9) 


0.4944(13) 


1.5 


4 


0.454134(11) 


0.6677(12) 


0.4923(19) 


1.4 


5 


0.454157(14) 


0.6711(16) 


0.4866(26) 


0.18 


6 


0.454160(15) 


0.6717(19) 


0.4857(31) 


0.16 


7 


0.454161(17) 


0.6717(23) 


0.4856(38) 


0.19 


8 


0.454164(20) 


0.6723(29) 


0.4846(49) 


0.22 


9 


0.454165(24) 


0.6724(38) 


0.4844(66) 


0.29 


10 


0.454177(45) 


0.6750(90) 


0.4796(160) 


0.38 


11 


0.454192(65) 


0.6783(136) 


0.4734(250) 


0.65 
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TABLE VIII. Results of the energy density E, the specific heat Ch, the susceptibihty x the 
second moment correlation length ^2"^ obtained from simulations at the fixed coupling K = 0.45417 
near the final estimate of the critical coupling, r denotes the integrated autocorrelation time of the 
specified observable, given in units of the average number of clusters needed to cover the volume 



of 


(ho lailic(\ 


The ^ 


itaiistics is given 


iu Icrnis 


of A" uioasurouionls Uikcu v 


very No upc 


late sl('|")s. 


L 


N 


No 


E 


TE 


Ch 


X 






4 


100k 


10 


0.40440(44) 


2.0(1) 


6.561(27) 


19.095(34) 


1.84(5) 


2.3104(37) 


8 


95k 


20 


0.35585(20) 


2.4(1) 


8.890(39) 


77.80(15) 


1.97(5) 


4.6852(65) 


16 


100k 


40 


0.338945(7) 


2.6(1) 


10.757(66) 


309.95(60) 


1.96(7) 


9.447(15) 


32 


83k 


80 


0.332815(3) 


3.1(1) 


12.520(73) 


1216.0(2.7) 


2.11(5) 


18.922(38) 


64 


72k 


160 


0.330628(2) 


3.7(1) 


14.35(11) 


4732(12) 


2.32(7) 


37.793(77) 



TABLE IX. Estimates for Kc{L) obtained via Binder's cumulant crossing technique of the 
reweighted fourth-order cumulants Ul and U1/2 nearest neighbour observable AW. Li — L2 
gives the pair of linear lattice sizes which determine the intersection point. 





Li — L2 


Ul 




Ul/2 


ATAT 


4-8 


0.4565(4) 




0.4617(3) 


0.4378(4) 


8-16 


0.4544(2) 




0.45457(8) 


0.4517(1) 


16-32 


0.45423(5) 




0.45424(4) 


0.45393(4) 


32-64 


0.45420(2) 




0.45421(2) 


0.45415(2) 



TABLE X. Results for the fourth-order cumulants Ul, Ul/2 the nearest neighbour ob- 



servable AW at = 0.45420 obtained with the reweighting technique from the simulations at 

K = 0.45417. The errors are obtained by a Jackknife-blocking procedure. 

L Ul Ul/2 NN 

4 0.59640(42) 0.56860(25) 0.70557(61) 

8 0.59134(42) 0.55270(31) 0.77439(45) 

16 0.58966(43) 0.55040(32) 0.79925(44) 

32 0.58907(50) 0.54974(37) 0.80640(48) 

64 0.58909(44) 0.54925(33) 0.80963(49) 
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TABLE XI. Estimates for the static critical exponent v obtained using eg. (|32|) , where A is 
replaced by the fourth-order cumulants C//, and ^"^^ the nearest neighbour observable AW 

with the critical coupling is set to = 0.45420, the final estimate of the critical coupling. 



lattice 




u 






L\ — L2 


Ul 




Ul/2 


AW 


4-8 


0.6496(93) 




0.5807(50) 


0.8443(86) 


8-16 


0.6799(111) 




0.6576(74) 


0.7519(76) 


16 - 32 


0.6649(126) 




0.6694(84) 


0.6977(68) 


32 - 64 


0.6584(154) 




0.6565(103) 


0.6779(83) 



TABLE XIL Estimates for the ratio of the static critical exponents j/u obtained using eq.(|3^. 
The first column gives the results of the ratio obtained from the susceptibility of the full lattice 
while the second column is obtained from the susceptibility of the subblocks. Li — L2 gives the 
pair of lattices, which is used to calculate the ratio of the exponents. 



lattice 




7/1/ 




Li — L2 


full lattice 




subblocks 


4-8 
8-16 
16 - 32 
32 - 64 


2.027(4) 
1.996(4) 
1.978(4) 
1.979(5) 




1.898(2) 
1.954(3) 
1.966(3) 
1.974(4) 


TABLE XIIL Expectation values of the specific heat Ch, the second moment correlation length 
^2"d, and the product of the helicity modulus T times the linear size of the system that are 
reweighted to the final estimate of the critical coupling Kc = 0.45420. The errors are calulated 
from a Jackknife analysis. 


L Ch 






T-L 


4 6.561(27) 
8 8.877(33) 
16 10.704(72) 
32 12.564(63) 
64 14.406(102) 




2.3112(34) 
4.6856(68) 
9.4639(142) 

19.003(34) 

38.25(69) 


1.090(2) 
1.091(4) 
1.12(1) 
1.13(2) 
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TABLE XIV. Comparison of critical properties determined from various methods. The results 
given for the simulations of the model in the high temperature phase are obtained from fits ac- 
cording to the simple scaling ansatz with the coupling K as parameter. For 7 and v from Ref. |25] 
we took the estimates of the fee lattice, since the errors are smaller than those obtained from the 
sc lattice. 



Met no a 


Kei. 




7 


V 


V 


T J 


Phenomenological RG 


this work 


U.4o42U(zj 




O.DD2(7j 


0.024(6) 


n con/ o^ 

0.589(2) 


TT" „,1 J- X. A /T/^ 

High temperature ML/ 


this work 


U.45417U(7j 


l.oz4(lj 








Phenomenological RG 


drill 


U.4o4z(ij 




0.670(zj 


,-^.n no 


0.586(1) 


High temperature MC 


§ 


0.45408(8) 


1.316(5) 








Phenomenological RG 


i 






^ 0.67 




0.590(5) 


High temperature MC 


i 


0.45421(8) 


1.327(8) 








e-expansion 


i 




1.315(7) 


0.671(5) 


0.040(3) 




e-expansion 


11 










0.552 


Resummed perturbation series 







1.3160(25) 


0.669(2) 


0.033(4) 




High temperature series 




0.4539(12) 


1.323(15) 


0.670(7) 


0.028(5) 




High temperature series 


m 


0.45414(7) 


1.325 


0.673 


0.030 




High temperature series 




0.45406(5) 


1.315(9) 


0.68(1) 






Experiment "^He 








0.6705(6) 
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